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Abstract 

Using semiclassical periodic orbit theory for a chaotic system in a 
weak magnetic field, we obtain the form factor predicted by Pandey 
and Mehta's two matrix model up to the third order. The third order 
contribution has a peculiar term which exists only in the intermediate 
crossover domain between the GOE (Gaussian Orthogonal Ensemble) 
and the GUE (Gaussian Unitary Ensemble) universality classes. The 
exact expression is obtained by taking account of the contribution 
from encounter regions where orbit loops are connected. 
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1 Introduction 



Universal level statistics of classically chaotic systems in agreement with the 
RMT (Random Matrix Theory) prediction has been intensively studied since 
it was conjectured two decades agojTj. 

Berry [2] first evaluated the first order term of the spectral form factor 
using the semiclassical periodic orbit theory 0. His method is called the 
diagonal approximation. It is applicable to both of the GOE (Gaussian 
Orthogonal Ensemble) and GUE (Gaussian Unitary Ensemble) universality 
classes. In the absence of time reversal invariance, the first order term is 
evaluated as a sum over the periodic orbits and agrees with the prediction 
derived from the GUE of random matrices. If the system is time reversal 
invariant, we additionally need to sum up over the pairs of the periodic 
orbits mutually time reverse. The result is also in agreement with the RMT 
prediction derived from the GOE. 

Recently Sieber and Richter (SR) showed that a family of orbit pairs 
with an encounter gave the first correction to the diagonal approximation 
jl]. Their idea clarifies what corresponds to the order of each term in the 
context of the periodic orbit theory. Heusler et al. |H] further generalized 
SR's idea by using a combinatorial analysis and Miiller et al. [HJ Ej obtained 
the full expansion which yields the RMT form factor. 

In this paper, we investigate the form factor of chaotic systems in a weak 
magnetic field, extending the method of Ref. 0. Our aim is to reproduce the 
universal level statistics in the crossover domain between the GOE and GUE 
universality classes. The corresponding RMT formula is derived from Pandey 
and Mehta's two matrix modeKS]. This formula interpolates the form factors 
of GOE and GUE, which correspond to chaotic systems without and with a 
magnetic field, respectively. Using a quantum graph, the authors previously 
obtained the second order term by employing Sieber and Richter's periodic 
orbit pairs and showed that it agreed with the RMT prediction [H]. Here we 
consider the third order term of the form factor and clarify a mechanism to 
yield the RMT expression. We show that the RMT expression is obtained by 
taking account of the interference of the gauge potential in the encounters of 
periodic orbits. 

Let us explicitly present the RMT prediction of the form factor. It is 
derived from Pandey and Mehta's two matrix model, as explained injD]. For 
small t(0 < t < |), the RMT form factor K^m{ t ) with a parameter /i is 
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written as 



K RM (r) = r + lf 1 dk-^—e-^ 



2 A-2r k + 2r 
= r + re - " - 2r 2 e^ + r 3 ^e" M + 2e~^ + • • • . (1) 

The GOE and GUE limits correspond to // — ► and /i — > oo, respec- 
tively. Up to the second order of r, each term is a monotonic function of 
the parameter \x. However, the third order term in (|TJ) includes a peculiar 
nonmonotonic term r 3i |-e _M , which vanishes in both of the GOE and GUE 
limits. Such terms generally appear in higher order contributions. We can 
thus expect that a semiclassical analysis of the third order term will offer a 
basis to reproduce the full expansion of Xrm(t)- 



2 Chaotic System in a Magnetic Field and 
Periodic Orbit Theory 

We consider a quantum system in a magnetic field and assume that the 
corresponding classical dynamics is chaotic. In the semiclassical limit H — > 0, 
the form factor of such a system can be formally written in terms of the 
classical periodic orbits [3 ; 

K po (r) = ^ Al A;^ +d ^e-« s y + WH (r - ^r^) ) , (2) 

where 5 7 , T 7 and A 7 are the action without gauge potential (including the 
Maslov phase), period and dimensionless stability amplitude of the periodic 
orbit 7. The phase 9 7 is an effect of the gauge potential in the action along 
the orbit 7, i.e., # 7 = / 7 <ix 7 ■ a(x 7 ) = f dtcbcy/dt ■ a(x 7 ), using the gauge 
potential a(x 7 ) at the position x 7 . The angular bracket stands for an average 
over the mean energy (E + E')/2 and over a time interval much smaller 
than the Heisenberg time Th- Since the classical dynamics is chaotic, each 
periodic orbit shows a diffusive behavior |1U[ ITT] : successive changes of the 
velocity can be regarded as independent events. Hence, if the time T elapsed 
on an orbit 7 is sufficiently large, we assume that the term g^(t) = d^/dt ■ 
a(x 7 ) can be replaced by the Gaussian white noise satisfying the correlation 
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((<7 7 (t)<7 7 (t'))) = 2D5(t — t'). Here a Gaussian average ((• • •)) is computed as 
a functional integral 



J Vg 1 exp 


AD J 


( T dt{ 9l {t)f 






J 


1 ° Vg 1 exp 


AD J 


( T dt{ 9l {t)f 







« %] )) 



Incorporating such a Gaussian average, we obtain a formula 

T-^iy ~~\~ Try 



(3) 



\ 7 ,7' 



2T H 



(4) 



We start from Berry's diagonal approximation, which yields the first order 
term [2]. This term is obtained from two types of periodic orbit pairs (7,7) 
and (7,7). Here a periodic orbit 7 is the time reverse of 7. For the pairs of 
identical orbits (7,7), the dependence on the gauge potential is canceled, so 
that the result does not depend on the magnetic field. The estimate for the 
form factor is given by the Hannay and Ozorio de Almeida (HOdA)'s sum 
ruler 



T 



T. 
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For the pairs of mutually time reverse orbits (7,7), on the other hand, the 
result depends on the magnetic field. It is estimated as 



J2 \A,\ 2 d(r-^) (( e ^-%)/*» = r((e 2i9 ^ h )) = re 



bT 



(6) 



where the last equality results from the Gaussian average (JSJ) . Here T is the 
period T(= tT#) and b is defined in terms of the diffusion constant D as 
b = AD/h 2 , which depends on the magnetic field. Putting the above results 
together, we reproduce the first order term of the form factor (THj as 



po 



r + re 



-bT 



(7) 



which reproduces the first two terms of the RMT prediction (fljl. The RMT 
parameter \i is identified with bT. 
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We next consider the second order term. As shown by Sieber and Richter 
(SR) JJ], a family of orbit pairs with one encounter plays a crucial role P ITU]. 
We can symbolically write SR pairs as 7 = EL\EL 2 and 7' = EL\EL 2 (see 
Figure 1). Here E and its almost time reverse E are the orbit segments in the 
encounter where two loops are connected. The loops are denoted as Li and 
L 2 , respectively, and L\ is the time reverse of L\. Let us see the behavior 
of the periodic orbits 7 and 7' in the encounter region. The directions of 

encounter orbits are depicted as Z^L for 7 and X for 7'. For a system 
with two degrees of freedom, it is convenient to use the Poincare section V 
in the phase space JSJ. The orbit 7 pierces through V at two phase space 
points denoted by vectors x a and x;,. Let us express Tx b — x a (T is the 
time reversal operator) in terms of the pairwise normalized vectors e s and e u 
as Txi, — x a = ue u + se s . To demand that the orbits in the encounter are 
mutually close, we introduce a bound c and assume that \u\, \s\ < c. Then 
the partner orbit 7' (or the time reverse of 7') pierces at x p a = x a + ue u and 
at x^ = T (x a + se s ). We can consequently estimate the duration t enc of the 
encounter and the action difference AS for the two orbits as 

1 c 2 

t C nc = — In - — r and AS" = us, (8) 
A \us\ 

where A is the Liapunov exponent. 

We need to estimate the number of encounters in one periodic orbit of a 
period T. This can be computed as 

re rT-2t cnc J 1 

/ duds / dT\—. — , (9) 

J-c Jo Nt cnc \l 

where f2~ 1 (= 2 ^ T ) is the ergodic return probability per a unit action. The 

factor T J T ~ 2tcnc dTi indicates that one of the two piercings occurs in the 
time interval [0, T] and that the time T\ elapsed on the loop L\ or L\ lies in 
[0, T — 2t enc \. The combinatorial factor J\f = 2 is necessary to avoid a double 
counting due to the equivalence of E and E in 7. 

Let us now evaluate the contribution of the gauge potential to the second 
order term. The Gaussian average yields a factor e~ bTl as a contribution 
from the the pair of the loops L\ and L\. In the region of the encounter, the 
gauge potential destructively interferes with itself and yields no contribution. 
That is, since mutually almost time reversed orbits are combined in the 
encounter, a product of the phase factors e t9 ~> and e~ ld ~<' gives a factor 1. 
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Incorporating the factor 2 due to the time reversal of the orbit 7', we obtain 
the second order term as 

-> -2rV bT , (10) 

again in agreement with Here we calculated only the term independent 
of £ enc , as the other terms vanish due to appearances of extra factors h or 
rapid oscillations [Zj. 



3 Third Order Term 

We now calculate the third order contribution, following the recipe explained 
in previous section. We will find that the calculation is somewhat different 
from the second order term because the interference of the gauge potential in 
the encounter regions gives a significant contribution. For the periodic orbit 
pairs which contribute to the third order term, Heusler et al.jH] introduced 
several notations (aas, api, ppi, ac and pc). In order to take account of the 
time reversal of the orbit 7', we modify their notations by putting suffixes 
(-a and -b) as drawn in Figures 2 and 3. 

Let us first consider a general orbit pair a with L loops and V encounters. 
Extending the result Q for the second order term, we can compute the 
number of encounters in one periodic orbit of a period T as 

rT—t a rT—t a —Ti rT—t a —Ti—T2 

duds / d7\ / dT 2 --- dT L _! Q a , (11) 



Jo Jo 



where 



T 

® a = Af U v t W-v' ^ 

- /v a llv=l b enc,v " 

The integration measures are du = Ilj=i duj and ds = Y[f=i dsj in terms 
of appropriate phase space coordinates (v,j, Sj). The times elapsed on the L 
loops are depicted as Ti,T 2 , ■ • ■ , T^. The duration of the f-th encounter is 
written as t enc ,v The total duration of the encounters in one periodic orbit 
is t a = J2v=i Wenc,« {h is the number of the orbit segments contained in the 
-y-th encounter). A combinatorial factor J\f a depends on the structure of the 
orbit pair. 
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We then consider the contribution of the gauge potential. The Gaus- 
sian average (J3J) on the loops yields a factor e~ bTl e~ bT2 ■ ■ -e~ bTK , where K 
is the number of the pairs of mutually time reversed loops. In Figures, the 
encounter orbits are schematically drawn along with the passing directions 
through the encounter regions. For the f-th encounter, we introduce the 
number n enc v of the orbit segment pairs which cause the significant gauge 
potential contribution. Let us fix an arbitrary direction (+) of the orbit 
passing and call the opposite direction (— ). Then the number n enCj „ can be 
computed as 



n, 



ciic,i; 



2 



{# {+) (7) " # H (7)} - {# {+) (7') " # ( - ) (7')}| • (13) 



Here # (+) (7) and #<-%) mean the numbers of passings of 7 through the 
encounter in (+) and (— ) directions, respectively. We again compute the 
Gaussian average © to evaluate the effect of the gauge potential. The contri- 
bution from an encounter is then estimated as ^ e 2mcntvu6,7//?l )) = e~ bricnc - vtcnc ' v . 

Putting the above results together, we find the contribution to the form 
factor from the orbit pair a as 

K p0) a{T) = t f duds 

x ' 



T—t a rT-t a -T\ rT-t a —Ti~T'2 T L _ 2 

dTi / dT 2 ■ ■ ■ / dT L . x Q a R a e ll ^ slh 



Jo jo 



(14) 



with 



R a = e -Kri+r2 + ---+T^) e -6(n2 nc>1 t en cA+^ ft c,2*e Il o,2+---+n^ CiV ,te n c,v)_ 



This contributes to the terms of order L — V + 1. Here r appears due to the 
HOdA sum rule and the action difference is estimated as AS = u j s j- 
Let us go back to the calculation of the third order term with L—V+l = 3. 
From Figures 2 and 3, we can readily see the followings. For a = aas,api and 
ppi, L = 4, V = 2 and t a = 2 (t en c,i + t cnCi 2)- On the other hand, for a = ac 
and pc, L — 3, V — 1 and t a = 3t enc . Thus we can calculate R a as 

#aas-a = #aas-b = e^ 1 ^, (16) 

D D „ — fr(7l+T 2 + tenc,l+£enc,2) (M\ 

Jl -api— a Jl - a pi— b c , \ LI J 

-Rppi-a = 1, (18) 
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&(T-tp pi )-4&(tenc,l+W,a) 

) 

&(Tl+T 2 ) 

) 

6(Tl+t cnc ) 

? 

-b(T-t pc )~9bt cnc 

Moreover the combinatorial factor Af a is evaluated in [Hj as 
(A4a S ,Ar api ,Ar ppi ,A4c,AT pc ) = (2,2,4,1,3). 



-Rppi— b 


= e 


-^ac— a 


= e~ 


Rac—b 


= e~ 


Rpc—a, 


= 1, 


Rpc—b 


= e" 



(19) 
(20) 
(21) 
(22) 
(23) 



(24) 



We substitute these formulas into (jHj) and again extract the terms inde- 
pendent of the encounter durations t e nc,i, 4nc,2 and t enc . After a straightfor- 
ward calculation, we obtain contributions from the orbit pairs 



-^"po,aas— a (T 



-^po,api— a (t 



po,ppi— a 



(r 



po.ppi 



-b [T 



-Apo.ac— a (T 



po,ac 



-b IT 



-Apo,pc— a(^~ 



K- 



po,pc- 



-b T 



po,aas- 



|)( s,api— b (t) 



A', 

r 3 



b(r) = 2r d e 

^3 



3-bT 



e -bT _ j ! + e -feT 



6T 



+ 6T + 1 



-3r^e 



3„-&T 



r 



— T 



2 e 



-6T 



6T 



-3 -&T 



(25) 
(26) 
(27) 
(28) 
(29) 

(30) 

(31) 
(32) 



-r-e " (1 + &T) . 

Summing up these values, we arrive at the third order term of the form factor 

'{bTf 



K${t) = Y j K po , c 



T = T~ 
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_ e -bT + 2e -bT 



(33) 



which agrees with the RMT prediction One might think that K poap [ 
and -Kpo,ac should disappear in the limit b —>■ oo, because a pair of mutually 
time reversed loops exists. However, as these contributions come from the 
limiting situations that the times elapsed on those loops are short, there is 
no inconsistency. 
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In summary, we investigate the level statistics of a classically chaotic 
system in the crossover domain between GUE and GOE universality classes. 
Summing up the contributions from the relevant periodic orbit pairs of a 
chaotic system, we find an agreement between the semiclassical form factor 
and the RMT result up to the third order. We expect that this agreement 
holds up to any arbitrary order. The periodic orbit pairs consist of loops and 
encounters. The durations of the encounters are logarithmically divergent 
in the limit H — > but infinitesimally small compared to the orbit periods, 
which are of the order of the Heisenberg time Th (~ 0(l/h) for a system 
with two degrees of freedom) j^j. Nevertheless, the contribution from the 
encounter regions is crucial. This mechanism should also be a key factor 
in the derivation of higher order terms. It might moreover shed light on a 
similar counting problem |14j for a quantum graph. 

Acknowledgement 

One of the authors (T.N.) is grateful to Dr. Gregory Berkolaiko, Dr. Petr 
Braun, Dr. Sebastian Miiller and Dr. Holger Schanz for valuable discussions. 

References 

[1] O. Bohigas, M.J. Giannoni and C. Schmit, Phys. Rev. Lett. 52 
(1984) 1. 

[2] M.V. Berry, Proc. R. Soc. London A400 (1985) 229. 

[3] M. Gutzwiller, Chaos in Classical and Quantum Mechanics 
(Springer, 1990). 

[4] M. Sieber and K. Richter, Physica Scripta T90 (2001) 128. 

[5] S. Heusler, S. Miiller, P. Braun and F. Haake, J. Phys. A37 (2004) 
L31. 

[6] S. Miiller, S. Heusler, P. Braun, F. Haake and A. Altland, Phys. 
Rev. Lett. 93 (2004) 014103-1. 

[7] S. Miiller, S. Heusler, P. Braun, F. Haake and A. Altland, 
|nlin. CD/05030521 



9 



[8] A. Pandey and M.L. Mehta, Commun. Math. Phys. 87 (1983) 449. 

[9] T. Nagao and K. Saito, Phys. Lett. A311 (2003) 353. 

[10] M. Turek and K. Richter, J. Phys. A36 (2003) L455. 

[11] K. Richter, Semiclassical Theory of Mesoscopic Quantum Systems 
(Springer, 2000). 

[12] J.H. Hannay and A.M. Ozorio de Almeida, J. Phys. A17 (1984) 
3429. 

[13] O. Bohigas, M.-J. Giannoni, A.M. Ozorio de Almeida and C. 
Schmit, Nonlinearity 8 (1995) 203. 

[14] G. Berkolaiko, H. Schanz and R.S. Whitney, J. Phys. A36 (2003) 
8373. 



10 



U Li (Li) 




SR pair 

Figure 1: The Sieber-Richter (SR) pair of periodic orbits 
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Figure 3: The periodic orbit pairs with one encounter. 
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